ON THE MASS FORMULA AND WIGNER AND 
CURVATURE ENERGY TERMS 



G. ROYER 

Laboratoire Subatech, UMR : IN2P3/CNRS-Universite-Ecole des Mines, 
4 rue A. Kastler, BP 20722, 44307 Nantes Cedex 03, France, 
E-mail: royer@subatech.in2p3.fr 

(Received February 2, 2008) 



Abstract. The efficiency of different mass formuias derived from the iiquid drop modei 
including or not the curvature energy, the Wigner term and different powers of the relative 
neutron excess 7 has been determined by a least square fitting procedure to the experimental 
atomic masses assuming a constant Ro.charge/A^'^'' ratio. The Wigner term and the curvature 
energy can be used independently to improve the accuracy of the mass formula. The different 
fits lead to a surface energy coefficient of around 17-18 MeV, a relative sharp charge radius 
ro of 1.22-1.23 fm and a proton form-factor correction to the Coulomb energy of around 0.9 
MeV. 

Key words: Nuclear masses, liquid drop model, Wigner term, curvature energy, charge 
radius. 



1 INTRODUCTION 

The binding energies of exotic nuclei close to the proton and neutron drip 
lines or in the superheavy element region are still poorly known and the dif- 
ferent predictions do not agree completely. Therefore continuous efforts are 
still needed to determine the nuclear masses. Within a charged liquid drop ap- 
proach, semi-macroscopic models including a pairing energy have been firstly 
advanced to reproduce the experimental nuclear masses Ul The coeffi- 
cients of the Bethe-Weizsacker mass formula have been determined once again 
recently [3J. To reproduce the irregularity of the masses as functions of A 
and Z partly due to shell closings and proton and neutron number parity, 
macroscopic-microscopic approaches have been elaborated, mainly the finite- 
range liquid drop model and the finite-range droplet model [1]. The Thomas- 
Fermi statistical model with a well-chosen effective interaction [5l |6] has also 
allowed to reproduced accurately the nuclear masses. Microscopic Hartree- 
Fock self-consistent theories using mean-fields and Skyrme or Gogny forces 
and pairing correlations [3 IB] as well as relativistic mean field calculations [S] 
have also been developed to describe these nuclear masses. Finally, nuclear 
mass systematics using neural networks have been undertaken recently jlOj. 

The evolution of the nuclear binding energy with deformation and rota- 
tion governs the fission, fusion, cluster and a decay potential barriers and the 
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characteristics of the large deformed rotating states. One or two-body shape 
sequences have to be selected to simulate the exit or entrance channels in 
the macroscopic- microscopic models. Within a generalized liquid drop model 
and a quasi-molecular shape sequence the main features of these barriers have 
been reproduced using, firstly, four basic macroscopic terms : the volume, 
surface. Coulomb and nuclear proximity energy contributions and, secondly, 
shell and pairing energy terms to explain structure effects and improve quan- 
titatively the results [121 [HI [HI [EJ HH [JTj . 

The purpose of the present work is to determine the efficiency of differ- 
ent combinations of terms of the liquid drop model by a least square fitting 
procedure to the experimentally available atomic masses [18] and to study, 
particularly, the separated influence of the Wigner term, the curvature energy 
and different powers of the relative neutron excess / to improve the GLDM. 

2 NUCLEAR BINDING ENERGY 

The nuclear binding energy B„„cK^i2) which is the energy needed to separate 
all the nucleons forming a nucleus is linked to the nuclear mass M.n.m by 

BnuciiA, Z) = Zmp + NniN - M„.„(^, Z). (1) 

This quantity may be connected to the experimental atomic masses given in 
[18] since : 

M„.„(A, Z) = Ma.m{A, Z) - Znie + Be{Z). (2) 



The binding energy Be(Z) of all removed electrons is determined by [19] 

B,{Z) = a,iZ^-^^ + b,iZ''-^\ (3) 

with aei = 1.44381 x 10"^ MeV and h^i = 1.55468 x lO^^^ ^eV. 

The following expansion of the nuclear binding energy in powers of A~^/^ 
and / = {N — Z)/A has been considered : 

. \ 1 3 e^Z^ Z^ 
-Ofc (1 - /cfcj/| - kkj'^ - kksl^) As - aoA^ + fp— 

^ ' 5rg^3 A 

— W\I\ + Epair — Eshell — Econg- (4) 

The first term is the volume energy and corresponds to the saturated 
exchange force and infinite nuclear matter. It includes the asymmetry energy 
term of the Bethe-Weizsacker mass formula via the I^A term. The second 
term is the surface energy term. It takes into account the deficit of binding 
energy of the nucleons at the nuclear surface and corresponds to semi-infinite 
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nuclear matter. In the Bethe-Weizsacker mass formula the dependence of the 
surface energy on / is not considered. The third term, the curvature energy, is 
a correction to the surface energy appearing when the surface energy is viewed 
as a result of local properties of the surface and consequently depends on the 
mean local curvature. This term is taken into account in the Lublin-Strasbourg 
drop (LSD) model [20], the TF model [6] but not in the FRLDM [4|. The 
term appears when the surface term of the liquid drop model is extended to 
include higher order terms in A~^/^ and /. The fifth term is the Coulomb 
energy. It gives the decrease of binding energy due to the repulsion between 
the protons. In the Bethe-Weizsacker mass formula the proportionality to 
Z(Z-l) is prefered. For the charge radius the formula Ro^charge = tqA^I^ is 
assumed, a more sophisticated expression has been studied in ref. [2T]. The 
Z'^jA term is a proton form- factor correction to the Coulomb energy which 
takes into account the finite size of the protons. The term proportional to / 
is the Wigner energy [H [23] which appears in the counting of identical pairs 
in a nucleus. 

The pairing energy has been determined using 

Epair = — Op/A-^/^ for odd Z, odd N nuclei, 

Epair = for odd A, (5) 
Epair = o.p/A^^'^ for even Z, even N nuclei. 

The Op = 11 value has been adopted following different fits. More sophisticated 
expressions exist for this pairing energy [U [6]. The theoretical shell effects 
given by the Thomas- Fermi model {7^^ column of the table in [5] and [6]) 
have been retained since they reproduce correctly the mass decrements from 
fermium to Z = 112 [22]. They are calculated from the Strutinsky shell- 
correction method and given for the most stable nuclei in the appendix. The 
sign for the shell energy term comes from the adopted definition in [5] . It gives 
a contribution of 12.84 MeV to the binding energy for ^''^Pb for example. The 
congruence energy term is given by : 

Econg = -lOMeV exp (-4.2|/|) . (6) 

It represents an extra binding energy associated with the presence of congruent 
pairs in contrast to the pure Wigner expression simply proportional to / [6]. 

The masses of the 2027 nuclei verifying the two following conditions 
have been used to obtain the coefficients of the different expansions by a 
least square fitting procedure : N and Z higher than 7 and the one standard 
deviation uncertainty on the mass lower than 150 keV [18]. The root-mean- 
square deviation has been calculated using : 

^2 ^ S [Mth - ME^pf 
n 
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In Table [H the improvement of the experimental mass reproduction when 
additional contributions are added to the basic ^, AP, A^/^, A^/^P, Z^/A^/^ 
terms is clearly shown (each calculation corresponds to one numbered line). 
The curvature energy is not taken into account. The introduction of the 
pairing term and of the proton form factor is obviously needful. In contrast, 
the congruence energy term does not allow to lower a at least with the fixed 
coefficients adopted here (as in the LSD and TF models). When the coefficients 
before the exponential and the exponent are free the congruence energy tends 
to the Wigner term since the coefficient before the exponential diminishes 
while the exponent increases. The constant term seems unnecessary. The 
term is useful to improve the accuracy of the expansion and is more 
efficient that the A2/3j4 r^^ie Wi gner term plays the major role to 

decrease a. When the Wigner term is taken into account the introduction of 
the A^/^\I\, and terms are ineffective. Thus, the very satisfactory 

value of cj = 0.60 MeV can be reached [H |8l [20]. The introduction of the 
Wigner term in a liquid drop model has the main drawback that it leads to an 
important discontinuity at the transition between one and two-body shapes 
as in fission or fusion. Indeed, when a single system divides into two parts 
the Wigner term must be evaluated separately for the two fragments and the 
results added. Thus for the same value of |/| (symmetric fission or fusion) 
the Wigner term will jump at scission to 2 times its original value. The same 
problem exists for the Congruence energy term. 

In Table O the efficiency of the curvature energy term with different / 
dependences is examined, disregarding the Wigner contribution. The intro- 
duction of only one term in A^^^ is ineffective while the addition of A^^^P 
improves slightly the results. Supplementary terms in |/| to determine the 
volume, surface and curvature energies allow to reach (7=0.59 MeV. They are 
still more efficient than P terms. The curvature energy term has the advan- 
tage that it is continuous at the scission point at least in symmetric fission. It 
has the disadvantage that its value (and the sign) lacks of stability. 

A good convergency of the volume a„ and asymmetry volume k^ con- 
stants is observed respectively towards around 15.5 MeV and 1.8 — 1.9 . The 
variation of the surface coefficient is larger but a^ evolves around 17-18 MeV. 
Small values of the surface coefficient favors quasi-molecular or two-body 
shapes at the saddle-point of the potential barriers while large values of 
promote one-body elongated shapes. The value of the proton form factor 
correction tends to 0.92 MeV (this value has been retained also in the LSD 
model). 

The reduced charge radius iCQ^charge converges to 1.22-1.23 fm. This is in 
full agreement with the set of 799 ground state nuclear charge radii presented 
in ref. ^24j. In this compilation a value of 0.9542A^/^ for the rms charge radius 
is obtained, which leads to 1.23yl^/'^ for the effective sharp charge radius. In 
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Tabic 1: Coefficient values (in MeV or fm) as functions of the selected term 
sets including or not the congruence and pairing energies and corresponding 
root mean square deviation. The shell energy is taken into account. 



n 


av 








as 


1 


14.8504 




1.55448 




16.1059 


2 


15.7826 


- 


1.6165 


- 


21.017 


3 


15.5959 


- 


1.70507 


- 


17.1723 


4 


15.6184 


- 


1.70459 


- 


17.242 


5 


15.5233 


- 


1.71612 


- 


18.071 


6 


15.4285 


- 


1.71066 


- 


17.5713 


7 


15.5718 


- 


1.60965 


1.96584 


18.0047 


8 


15.5919 


-0.04582 


1.90751 


- 


17.8069 


9 


15.447 


_ 


1.84011 


_ 


17.3581 


10 


15.4647 


-0.0049 


1.8509 


- 


17.3894 


11 


15.5147 


- 


1.85055 


- 


17.6976 


12 


15.4607 


- 


1.83014 


0.14095 


17.3706 


n 


ksi 






ao 


W 


1 




0.93696 




_ 


_ 


2 


- 


1.13845 


- 


-18.253 


- 


3 


- 


0.98894 


- 


- 


- 


4 


- 


0.98033 


- 


- 


- 


5 


- 


1.40391 


- 


- 


- 


6 


- 


1.39267 


- 


1.8569 


- 


7 


- 


0.75244 


10.56741 


- 


- 


8 


-0.283 


2.39131 


- 


- 


- 


9 


- 


2.11347 


- 


- 


27.5488 


10 


-0.03409 


2.16959 


- 


- 


25.1137 


11 


- 


2.14699 


- 


-1.409 


29.0793 


12 


- 


2.01948 


1.10833 


- 


26.6425 


n 


Cong 


Pairing 




fp 


a 


1 


y 


y 


1.2434 


2.52888 


1.156 


2 


y 


y 


1.1595 


3.34645 


0.936 


3 


n 


n 


1.2272 




1.322 


4 


n 


y 


1.2244 




1.032 


5 


n 


y 


1.2066 


1.47705 


0.687 


6 


n 


y 


1.2153 


1.39388 


0.684 


7 


n 


y 


1.2042 


1.24476 


0.665 


8 


n 


y 


1.2049 


1.0255 


0.628 


9 


n 


y 


1.2252 


0.95419 


0.603 


10 


n 


y 


1.2234 


0.93536 


0.603 


11 


n 


y 


1.2195 


0.98825 


0.601 


12 


n 


y 


1.2244 


0.91559 


0.602 



Table 2: Coefficient values (in MeV or fm) as functions of the selected term sets 
and corresponding root mean square deviation. The shell and pairing energies 
are taken into account but not the congruence energy and the Wigner term. 



n 


av 


k 

"'til 


k 

"'V2 


k 




1 


15.3416 


- 


1.70872 


- 


16.8356 


2 


15.3225 


- 


1.87616 


- 


16.9627 


3 


15.2954 


- 


2.39856 


-7.873 


16.8264 


4 


15.5668 


0.17993 


1.28391 


- 


18.5295 


n 


ksi 










1 




1.40978 




1.99142 




2 




2.94224 




1.2545 




3 




7.85795 


-73.9552 


0.709782 




4 


1.46203 


-2.1873 




-2.72953 


25.5927 


n 






ro 


fp 


a 


1 






1.2205 


1.33865 


0.68 


2 


-50.7382 




1.2281 


1.28911 


0.66 


3 


-395.7891 


4557.60 


1.2391 


1.03047 


0.61 


4 


-62.900 




1.2285 


0.91998 


0.589 



this adjustment to the nuclear masses the nuclear mass radius is not fitted. 
Root-mean-squared matter radii are given in ref. |25] for specific nuclei. 
For the Bethe-Weizsacker formula the fitting procedure leads to 

Bnud{A, Z) = 15.69A - 17.6037^2/^ (8) 

-0.71660^^^|j^ - 23.6745/2^ + E^air - Esheii 

with £7=1.35 MeV. That leads to ro=1.2057 fm and k^=1.5089. The non 
dependence of the surface energy term on the relative neutron excess / explains 
the a value. 

The Fig.[T]shows the dispersion between the theoretical and experimental 
masses within the last formula presented in Table [2] and given below : 



B = 15.5668 (l - 0.17993|/| - 1.283917^] A 
-18.5295 (l - 1.46203|/| + 2.1873/^) 
-f2.72953 f 1 - 25.59271/1 + 62.9/^) A^ 



3 e^Z^ 

■3 -+0.91998— +Epa^r-E,hea■ (9) 

5 1.2285^3 A 



6 



-1 




Figure 1: Difference (in MeV) between the tlieoretical and experimental 
masses for the 2027 nuclei as a function of the mass number. 



3 SUMMARY AND CONCLUSION 

The efficiency of different mass formulas derived from the liquid drop model 
and including or not a curvature energy term, the Wigner term and different 
powers of the relative neutron excess / has been determined by a least square 
fitting procedure to 2027 experimental atomic masses assuming a constant 
Ro^c/iarge/A^/^ ratio. The Wigner term and the curvature energy term can 
improve independently the accuracy of the mass formula. The very satisfactory 
value of cr = 0.59 MeV can be reached. The different fits lead to a volume 
energy coefficient of around 15.5 MeV, a surface energy coefficient of around 
17-18 MeV, a relative charge radius ro of 1.22-1.23 fm and a proton form-factor 
correction of around 0.9 MeV. The addition of a term in |/| in the volume, 
surface and curvature energy terms is more efficient than a term in I^. 
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Table 3: Theoretical shell energy (in MeV) extracted from [5] {7^^ column) at 
the ground state of nuclei for which the half-life is higher than 1 ky. 





1^0 


180 


I9p 




2iA^e 


22Are 


23Ara 


-0.45 


1.19 


1.3 


2.76 


2.81 


2.82 


2.19 


2.22 










2M/ 


28 5i 


29 5i 


305i 


1.64 


1.77 


0.64 


1.89 


0.79 


-0.26 


-0.25 


0.22 


31p 


32^ 


335 


345 


36^ 


35 c/ 


36 C7/ 


37C7Z 


0.22 


0.66 


0.86 


1.13 


1.22 


1.32 


1.48 


1.40 






40^r 


39^ 


40^ 


41^ 


40Ca 


4iC7a 


1.57 


1.64 


2.45 


1.79 


2.58 


2.58 


1.71 


2.49 






44Ca 


46C7o 


48Ca 


455c 


46ri 




2.49 


2.16 


1.7 


0.69 


-0.82 


2.43 


2.44 


1.95 








50 y 


51 y 


50 Cr 


52Cr 


53Cr 


1.44 


0.81 


-0.05 


0.54 


-0.31 


0.74 


-0.69 


-0.38 






55 Mn 


54pg 


56 Fe 


57pg 


58pg 


60pg 


0.74 


-1.11 


0.68 


-1.54 


0.05 


0.72 


1.22 


2.07 


59Co 


'^Ni 


'^Ni 




^^Ni 


62iVi 


64iVi 


63 Cn 


0.77 


-1.58 


-0.68 


-0.16 


0.59 


1.1 


1.63 


1.86 


^'Cu 




66 Zn 


67 Zn 


68 Zn 


™Zn 


69 Ga 


71 Ga 


2.33 


2.53 


2.89 


3.16 


2.99 


2.94 


3.79 


3.71 


™Ge 




^^Ge 


^^Ge 


76 Ge 




745e 


76 5e 


4.13 


4.08 


4.19 


3.82 


2.53 


4.08 


4.42 


4.08 






™5e 




82 5e 


79p^ 


8iSr 


80i^r 


4.06 


3.27 


2.87 


1.89 


0.38 


4.07 


2.28 


4.39 






83^^ 


^^Kr 


86 


85i?6 


87i?6 


86 5r 


3.77 


2.74 


1.65 


0.96 


-0.40 


1.13 


-0.35 


0.79 




88^^ 


89 y 


90 zr 


91 Zr 


92 Zr 


93 Zr 


94Zr 


0.05 


-0.97 


-1.19 


-1.63 


-0.47 


0.46 


1.53 


2.54 






93JV6 




92Mo 


93 Mo 


94Mo 


95 Mo 


3.49 


-0.57 


0.44 


1.51 


-2.12 


-1.07 


-0.12 


0.97 






98Mo 


ioomo 


s^Tc 


98rc 




'^Ru 


1.79 


2.46 


2.98 


3.62 


1.26 


2.05 


2.64 


-1.11 










i02i?n 


104i?^, 


103^/1 


i02pd 


0.57 


1.36 


2.00 


2.54 


2.98 


3.49 


2.44 


0.63 




105p^ 


106pj 




108p^ 


iiopd 


10M5 


109^^ 


1.83 


2.39 


2.80 


3.08 


3.34 


3.42 


2.20 


2.91 










112Cd 


ii3(:7d 


114Cd 


ii6Gd 


0.31 


1.35 


2.11 


2.42 


2.52 


2.61 


2.50 


2.26 


113 


115j^ 






ii^Sn 


ii65n 


ii75n 


ii85n 


1.83 


1.97 


0.37 


0.81 


1.03 


0.94 


0.96 


0.75 








i245n 


126 5n 


12156 


12356 


i20re 


0.70 


0.19 


-0.99 


-2.51 


-4.36 


0.76 


-0.16 


2.08 









l^^Te 


I2«re 


128pg 


i3Ure 


1-27 J 


1.55 


1.30 


0.66 


0.25 


-0.62 


-2.46 


-4.74 


0.35 


129 J 




i29Xe 


i30Xe 


i3iXe 


i32Xe 


134^6 


i36Xe 


-1.29 


1.01 


0.48 


-0.31 


-1.11 


-2.36 


-4.92 


-7.2 






1325a 


i34pa 


i35pa 


i36pa 


i37pa 


i38pa 


-1.28 


-3.90 


0.93 


-0.55 


-1.45 


-3.01 


-4.18 


-5.29 




i38La 


i39La 


i36Ce 


i38Ce 


i40Ce 


i42Ce 


141 p^ 


-2.22 


-3.37 


-4.50 


0.70 


-1.57 


-3.86 


-2.06 


-3.26 






144 iVd 


145 A^d 


146 iVd 


148 iVd 


150 A^d 


1445m 


-2.88 


-2.14 


-1.04 


0.10 


0.56 


0.85 


0.54 


-2.28 




i47Sm 


i48Sm 


1495m 


i™Sm 


152 Sm 


1545m 




-0.41 


0.67 


1.12 


1.22 


1.31 


0.90 


0.38 


1.39 






i52Gd 


i54Gd 


i55Gd 


i56Gd 


is^Cd 


158Gd 


1.02 


1.30 


1.59 


1.33 


1.05 


0.89 


0.62 


0.56 






154 


156p)y 


158^y 


160p,y 


161 


162 


0.21 


0.64 


1.63 


1.56 


1.24 


0.92 


0.64 


0.47 




164£,y 


165 


163i/o 


165//0 


162^^ 


164^^ 


166 p^ 


0.16 


-0.06 


-0.41 


0.46 


-0.12 


1.20 


0.70 


0.07 




168^^ 


170^^ 


i69rm 


168y^ 


170y^ 


171y^ 


172y^ 


-0.37 


-0.54 


-1.08 


-0.60 


0.32 


-0.34 


-0.76 


-0.94 




174y^ 


176^5 




176L^i 


174p-_^ 


176/// 


177p^ 


-1.32 


-1.30 


-1.74 


-1.23 


-1.62 


-0.38 


-0.90 


-1.33 




179^^ 


180^^ 


I82i7/ 


i«ira 


180 


182^ 


183 


-1.53 


-1.97 


-1.99 


-2.16 


-2.02 


-1.21 


-1.71 


-2.00 




186^ 


i85i?e 


i87pe 


i840s 


1860s 


is^Os 


1880s 


-2.02 


-2.38 


-2.19 


-2.48 


-1.61 


-1.88 


-2.16 


-2.08 




i90Os 


1920s 


191/r 


193 


IQOpi 


192 Pi 


194pi 


-2.43 


-2.47 


-3.50 


-2.54 


-3.62 


-0.97 


-2.01 


-3.31 


195 Pi 


196 Pi 


198 Pi 


19Mn 


196^^ 


198^^ 


199 p-^ 


200p^ 


-4.04 


-4.80 


-6.11 


-5.56 


-4.51 


-5.99 


-6.75 


-7.52 


201^^ 




204/7^ 




205 ji^ 


202p6 


204p6 


205 P6 


-8.37 


-9.11 


-10.69 


-9.97 


-11.58 


-8.22 


-10.02 


-11.00 


206p^ 


207p^ 


208 p^ 


208p^ 


209 p^ 


226 Pa 


229r/l 


230T/1 


-11.82 


-12.68 


-12.84 


-11.70 


-11.95 


-0.30 


-0.52 


-0.43 




231 Pa 


233 


234 


235^ 


236^ 


238^ 


236 ATp 


-0.60 


-0.79 


-1.27 


-1.23 


-1.46 


-1.30 


-1.27 


-1.85 


237^^ 


239p^ 


240pn 


242pu 


244pu 


243^^ 


245 Crn 


246Cm 


-1.74 


-2.12 


-1.95 


-1.99 


-2.08 


-2.44 


-3.05 


-2.96 




248Cm 


247pA; 












-3.17 


-3.00 


-3.46 
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